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SUMMARY 

A  computer  code  has  been  written  which  utilizes  ray-tracing  techniques  to 
predict  the  changes  in  position  and  geometry  of  a  laser  doppler  velocimeter 
probe  volume  resulting  from  refraction  effects.  The  code  predicts  the  posi¬ 
tion  change,  changes  in  beam  crossing  angle,  and  the  amount  of  uncrossing  that 
occur  when  the  beams  traverse  a  region  with  a  changed  index  of  refraction, 
such  as  a  glass  window.  The  code  calculates  the  changes  for  flat  plate,  cylin¬ 
der,  general  axisymmetric,  and  general  surface  windows  and  Is  currently  opera¬ 
tional  on  a  VAX  8600  computer  system. 


INTRODUCTION 

The  laser  anemometer  has  provided  the  fluid  dynamicist  with  a  powerful 
tool  for  nonintrusively  measuring  fluid  velocities.  One  of  the  more  common 
types  of  laser  anemometers,  the  laser  fringe  anemometer,  divides  a  single 
laser  beam  into  two  parallel  beams  and  then  focuses  them  to  a  point  in  space 
called  a  "probe  volume"  (PV)  where  the  fluid  velocity  is  measured. 

Many  applications  using  this  method  for  measuring  fluid  velocities  require 
the  observation  of  fluids  through  a  window.  Whenever  light  traverses  a  region 
of  a  different  refractive  index,  its  direction  of  travel  is  changed  in  a  man¬ 
ner  described  by  Snell's  law: 


N  sin  a  =  N'sin  a' 


(1) 


where 

N,N‘  indices  of  refraction  before  and  after  the  surface 

a, a'  angles  between  the  surface  normal  and  the  light  beam  on  either  side  of 
the  surface. 

With  the  exception  of  a  few  very  specific  window/optics  configurations, 
the  Implications  of  this  law  to  the  laser  probe  volume  when  a  window  is 
Inserted  In  the  laser  beam  optical  path  between  the  focusing  lens  and  the  probe 
volume  can  be  summarized  as  follows: 
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(1)  The  position  of  the  actual  probe  volume  will  change. 


(2)  The  beams  will  uncross,  i.e.,  no  longer  lie  in  the  same  plane  at  the 
probe  volume  location. 

(3)  For  nonflat  plate  windows,  the  crossing  angle  of  the  two  beams  will 
change. 

The  first  result  means  that  data  will  not  be  acquired  in  the  desired  loca¬ 
tion.  The  second  affects  the  data  acquisition  rate  by  distorting  the  probe 
volume,  which  only  exists  where  the  beams  cross.  The  last  impacts  the 
accuracy  of  the  measurements  taken  because  it  is  directly  related  to  the  calcu¬ 
lation  of  velocity  by  the  following  equation: 


V  =  Xf/<2  sin(e/2)>  (2) 


where 

X  light  wavelength 
f  fringe  crossing  frequency 
e  beam  crossing  angle 

In  the  past,  most  researchers  have  assumed  that  the  resulting  errors  are 
small  and  can  be  ignored.  The  few  attempts  to  assess  the  seriousness  of  win¬ 
dow  refraction  effects  (refs.  1  to  4)  have,  unfortunately,  been  concerned  with 
special  cases  such  as  flat  windows  and  cylinders  with  observations  in  a  single 
plane  normal  to  the  axis  of  symmetry.  This  approach  has  proved  reasonable  in 
the  past  since  windows  were  usually  simply  shaped  and  thin.  Further,  observa¬ 
tion  angles  were  generally  near  normal  to  the  window  surfaces. 

A  more  general  approach  developed  at  the  NASA  Lewis  Research  Center  to 
meet  the  needs  resulting  from  the  use  of  more  complex  window  shapes  and  pas¬ 
sage  geometries  is  reported  here. 

A  ray  tracing  technique  which  is  not  restricted  to  special  cases  is  used 
to  study  the  changes  in  probe  volume  geometry  and  position  due  to  refraction 
effects  caused  by  both  flat  and  general  smooth  windows.  This  technique  pre¬ 
dicts  the  laser  beam  paths  through  any  of  four  different  window  shapes: 

(1)  Flat  plate  windows  (constant  thickness) 

(2)  Cylindrical  windows  (constant  thickness) 

(3)  General  axisymmetric  windows  (constant  thickness) 

(A1*  General  smooth  windows  (variable  thickness) 

When  describing  the  window  surfaces,  the  technique  calculates  the  flat 
plate  and  cylindrical  window  cases  analytically.  The  general  axisymmetric 
case  is  treated  analytically  in  the  R-9  plane  and  uses  a  cubic  spline  fit  in 
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the  R-Z  plane  where  the  Z  axis  represents  the  axis  of  window  symmetry. 

The  smooth  general  window  case  uses  the  cubic  spline  fit  in  both  the  X-Y 
and  X-Z  planes. 

The  technique  yields  the  new  probe  volume  position,  the  new  crossing 
angle,  and  the  amount  of  beam  uncrossing  (specified  as  both  the  number  of 
probe  volume  fringes  and  the  absolute  minimum  distance  between  beams). 

Two  separate  computer  programs  are  used.  The  first  program  generates  an 
input  data  set  for  use  by  the  second  program.  The  second  calculates  the 
changes  in  the  geometry  of  the  probe  volume.  The  use  of  a  separate  program 
allows  the  creation  of  permanent  input  datasets. 

The  following  parameters  may  be  varied  to  describe  the  case  to  be 
evaluated: 

(1)  Indices  of  refraction  on  both  sides  of  the  window  and  the  window 
itself 

(2)  The  window  shape 

(3)  The  assumed  position  of  the  probe  volume  and  the  actual  position  of 
the  focusing  lens  relative  to  the  window 

(4)  The  orientation  of  the  plane  which  contains  the  laser  beams 

(5)  The  beam  crossing  angle 

(6)  The  laser  beam  wavelength 

This  report  contains  a  description  of  the  problem  geometry,  the  method  of 
analysis,  complete  information  for  using  these  programs,  example  input  and  out 
put,  and  several  example  cases  for  general  axi symmetric  and  smooth  totally  gen 
eral  surfaces. 

These  FORTRAN  77  programs  are  currently  operational  on  a  DEC  VAX8600  com¬ 
puter.  These  Interactive  programs  require  a  typical  core  memory  of  86  pages 
and  have  typical  CPU  run  times  of  less  than  1  sec,  depending  upon  window 
geometry. 


GEOMETRY 

A  basic  understanding  of  the  window  and  optics  geometry  and  the  variables 
used  to  describe  them  is  necessary  to  correctly  use  the  technique  explained  in 
this  report.  Therefore,  we  will  take  a  more  detailed  look  at  the  geometries 
of  the  surfaces  involved.  Figure  1  provides  information  on  the  general  orien¬ 
tation  of  the  probe  volume  and  lens  for  all  windows.  While  most  of  the  infor¬ 
mation  presented  in  this  figure  is  self  explanatory,  some  points  should  be 
clarified. 

The  X',  Y ' ,  Z'  coordinate  system  shown  in  figure  1  has  its  origin  at 
the  probe  volume.  The  X-Y-Z  coordinate  system,  not  shown  in  figure  1,  is 
oriented  such  that  the  Z  axis  is  along  the  axis  of  symmetry  of  the  window 
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being  described,  should  this  window  be  symmetric.  The  axes  of  the  X‘,  V, 

V  coordinate  system  are  parallel  to  their  respective  X,  Y,  Z  axes.  The 
origins  of  the  two  coordinate  systems  are  not  necessarily  coincident. 

The  X  axis  should  be  oriented  roughly  so  that  the  optical  system  bisect¬ 
ing  line  parallels  the  X-Z  plane. 

The  bisecting  line  is  defined  as  a  line  equidistant  from  both  laser  beams 
and  lying  in  the  plane  defined  by  them.  It  connects  the  lens  center  and  the 
original  probe  volume. 

The  original  probe  volume  is  the  probe  volume  position  that  would  exist 
if  there  were  no  window  in  the  optical  path  between  the  focusing  lens  and  the 
probe  volume. 

The  beam  crossing  angle,  e,  is,  for  all  cases,  the  angle  defined  as: 
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tan  e  =  | A]  x  A2  | /A]  •  A2 


(3) 


where 


beam  crossing  angle 


Ai ,  A2  incoming  laser  beam  vectors 

The  angle  4»,  not  shown  in  figure  I,  is  the  angle  between  the  surface  nor¬ 
mal  and  the  bisecting  line  for  the  flat  plate  window. 

For  all  cases,  the  beam  plane  orientation  angle,  6,  is  defined  as  the 
angle  measured  counterclockwise  in  the  lens  plane,  made  by  the  line  connecting 
the  two  beams  where  they  cross  the  lens  plane  and  a  line,  Z",  where  Z"  is  a 
line  lying  in  a  plane  formed  by  the  Z'  axis  and  the  beam  bisecting  line 
(fig.  1). 

Figure  2  shows  the  geometry  for  the  cylindrical  window  case  and  additional 
information  for  describing  the  lens  and  original  probe  volume  position.  Notice 
that  the  axis  of  symmetry  of  the  window  is  the  Z-axis. 

For  all  window  cases,  the  lens  center  and  probe  volume  position  are 
described  in  the  R,  Z,  9  coordinate  system.  Neither  the  probe  volume  nor 
the  lens  center  axial  position  are  restricted  to  a  Z  equal  zero  requirement. 

Figure  3  shows  the  general  axisymmetric  coordinate  system.  The  window  is 
described  as  a  set  of  R-Z  points  that  specify  the  window  outer  surface.  The 
window  inner  surface  is  specified  using  a  normal  thickness.  For  this  window, 
the  thickness  must  be  constant. 

Figure  4  indicates  the  method  in  which  the  general  smooth  window  is  speci¬ 
fied.  Here  both  surfaces  are  specified  as  an  array  of  X,  Y,  Z  points  in 
space.  There  are  no  restrictions  on  window  thickness  in  this  case.  The  num¬ 
bers  in  parentheses  indicate  the  array  locations  of  the  points. 
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The  specific  restrictions  to  the  geometries  described  above  are: 

(1)  Surfaces  must  be  smooth  and  smoothly  changing. 

(2)  For  all  general  axisymmetric  surfaces,  surface  slopes  in  the  Z 
direction,  DX/OZ,  must  remain  less  that  45°  and  DY/DZ  must  remain  less  than 
45°  for  all  general  smooth  surfaces. 

(3)  The  Z  axis  is  defined  as  the  center  of  rotation  for  the  cylindri¬ 
cal  and  general  axisymmetric  surface. 

(4)  R1 ,  the  lens  radial  position,  should  be  greater  than  Rpv,  the  probe 
volume  radial  position,  and  greater  than  Ro,  the  window  outer  radius  and 

Rpv  should  be  less  than  Ri ,  the  window  inner  radius  (fig.  2). 


METHOD  OF  ANALYSIS 


The  code  uses  a  simple  ray  tracing  technique  for  the  necessary  calcula¬ 
tions.  To  accomplish  this,  the  initial  ray  direction  cosines  and  the  points 
in  space  where  the  beams  cross  the  lens  plane  are  calculated. 

Then,  from  the  window  surface  Information,  the  point  of  interception  of 
each  beam  with  the  upper  or  outer  window  surface  and  the  surface  window  normal 
at  that  point  are  calculated.  Snell's  law  is  applied  to  each  beam  as  it 
crosses  the  surface.  The  same  procedure  is  then  applied  to  find  the  intercep¬ 
tion  of  the  beam  and  the  lower  or  inner  window  surface. 


With  the  position  and  the  direction  cosines  of  the  two  beams  at  the  inside 
of  the  window,  the  calculation  of  minimum  crossing  distance,  crossing  angle, 
and  the  point  of  nearest  approach  become  relatively  straightforward  geometry 
problems . 

Of  importance  is  how  the  beam/window  interception  point  and  the  surface 
normals  at  the  Interception  point  are  calculated. 

The  following  section  discusses  the  methods  used  to  determine  the  point 
of  beam/window  Interception  and  the  window  outward  surface  normal  at  the 
beam-window  interception  point  for  each  of  the  four  window  types  this  computer 
code  solves. 

At  the  point  of  intersection  of  the  laser  beams  and  surfaces,  the  laser 
beams  and  the  surface  normals  can  be  described  using  direction  cosines  in  the 
fol lowi ng  way: 


V  =  (cos  a) i  +  (cos  |3)j  +  (cos  y)k 


i 

(4) 


where 

V  incoming  laser  beam  vector 


a,  (3,  y  direction  angles  between  the  vector  and  the 


X,  Y,  Z  axes 


".>_i 


for  each  incoming  laser  beam  and 

N  =  (cos  u)i  +  (cos  Oj  +  (cos  o)k  (5) 

where 

N  outward  surface  normal 

v.C.o  direction  angles  between  the  vector  and  the  X,  Y,  Z  axes 

for  each  outward  surface  normal  (figs.  5  and  6),  where  both  are  defined  as  unit 
vectors . 

The  angle  between  these  two  vectors  is  simply  the  arc  cosine  of  the  dot 
product. 
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COS-14>  =  V  •  N  =  COS  a  COS  u  +  COS  [$  COS  C  +  COS  Y  COS  a 


(6) 


Using  a  simple  ray  tracing  analysis  and  Snell's  law;  the  direction  of  the 
ray  after  crossing  the  surface  can  be  determined  by  (ref.  5): 

V'  =  (n/n 1 >V  +  (b/n ' )N  (7) 


where  n  Is  the  index  of  refraction  on  the  incoming  side  of  the  surface,  n1 
is  the  index  of  refraction  after  crossing  the  interface,  N  is  defined  by 
Eq.  5,  and  b  is; 


b  =  n 1  cos  <j>  -  n  cos  <j>* 


(8) 


and  4>'  is  given  by 


COS  <t> '  =  (n'2  -  n2  +  n2  cos2  a)/n' 


(9) 


Using  these  equations,  a  method  for  describing  the  window  surfaces  in 
space,  and  the  window  surface  normals;  the  path  of  each  beam  can  be  traced. 

The  window  surfaces  and  their  surface  normals  are  described  differently 
for  each  of  the  four  types  of  surfaces. 

The  flat  window  surfaces  are  described  using  the  equation  of  a  flat 
surface : 


(X  -  XD/COS  V  =  (Y  -  YD/COS  c  =  (Z  -  ZD/COS  a  (10) 
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where 

X, Y,Z  surface  coordinates 

XI. Y1.Z1  a  specified  surface  point 
v.C.o  surface  normal  direction  cosines 


A  flat  surface  can  be  specified  by  any  three  points  not  in  a  straight 
line.  The  cross  product  of  any  two  unit  vectors  connecting  these  points  pro¬ 
vides  the  surface  normal  and  its  direction  cosines.  Since  the  surface  is 
flat,  the  surface  normal  has  a  constant  direction. 


The  cylindrical  and  general  axl symmetri c  windows  are  described  using  the 
equation: 


R2  =  X2  +  Y2 


(11) 


where 

R  constant  radius  for  a  cylinder  and 

R  f(z)  for  the  general  axisymmetric  case  and  f(z)  is  specified  by  a  series 

of  R-Z  points 

X, Y  surface  coordinates 

The  necessary  direction  cosines  for  a  cylinder  are  also  relatively 
straight  forward  calculations.  The  slope  of  the  circle  at  any  point  can  be 
obtained  by  taking  the  derivative  of  X  with  respect  to  Y.  This  informa¬ 
tion,  coupled  with  the  fact  that  the  direction  cosine  with  respect  to  the  axis 
of  symmetry  is  zero  provides  all  the  necessary  information. 

To  approximate  the  unknown  R-Z  relationship  in  the  general  axisymmetric 
case,  the  program  makes  use  of  a  cubic  spline  fit  subroutine  described  in  ref¬ 
erence  6.  It  solves  a  tridiagonal  matrix  given  in  reference  7.  The  inputs 
are  a  series  of  R-Z  values.  It  should  be  noted  that  the  spline  subroutine 
requires  that  the  end  point  condition  of  the  second  derivative  at  either  curve 
end  point  be  specified.  The  input  default  gives  an  end  point  second  deriva¬ 
tive  of  one  half  the  second  derivative  of  the  adjacent  spline  point.  Care 
must  be  taken  to  insure  the  input  points  accurately  describe  the  surface,  that 
the  surface  is  smooth  and  that  the  end  point  second  derivatives  are  not  greatly 
in  error.  Inaccuracies  occur  in  the  solution  when  the  surface  described  has 
slopes  approaching  90°  and  surfaces  of  this  type  should  be  calculated  using 
the  general  surface  routine  with  an  appropriate  coordinate  rotation  for  the 
geometry. 

Output  from  this  cubic  spline  fit  at  the  desired  Z  location  includes 
the  value  of  the  radius  and  both  the  first  and  second  derivatives  of  radius 
with  respect  to  axial  position.  To  solve  for  the  intersection  between  the  win¬ 
dow  surface  and  incoming  laser  beam,  the  incoming  beam  is  first  described  as  a 
function  of  radius  and  axial  position.  An  iterative  procedure  is  then  used  to 
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find  the  R-Z  interception  point  of  the  incoming  ray  and  the  surface.  The  de¬ 
sired  solution  accuracy  for  the  iterative  procedure  is  an  input  value. 

At  the  calculated  Intersection  point,  the  slopes  of  the  general  axisymmet- 
ric  surface  are  known: 

dR/dZ  and  dX/dY 

However,  we  are  interested  in  the  surface  normal  direction  cosines  (fig.  5): 
cos  v,  cos  C,  and  cos  a.  The  last  is  directly  available  since: 


tan  a  =  -1 /(dR/dZ)  (12) 


Now,  recalling  that: 

1  =  (cos2v)  +  (cos2C)  +  (cos 2o)  (13) 


and  that,  since  dX/dY  is  known 


tan  5  =  dX/dY  (14) 


Now, 


cos  v  =  b/s  (15) 


and 


s  =  (a2  +  b2  +  d2)1/2  (16) 


and  realizing  that 


a  =  b  tan  £  and  d  =  b  tan  <  (17) 


where  a,  b,  c,  and  d  are  distances  defined  in  figure  5(b)  we  can  write 

cos  v  =  b/(b2  +  b2tan2^  +  b2tan  k)1/2  (18) 


or 
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cos  v  =  1/(1  +  tan2£  +  tan2*)!^ 


(19) 


and 


cos  t;  =  (1  -  cos 2o  -  cos2v)l/2  (20) 


The  general  surface  solution  follows  the  same  logic;  the  difference  being 
that  no  analytical  expressions  exist  to  describe  the  surface.  The  surfaces 
are  described  by  an  array  of  X,  Y,  Z  coordinates.  These  coordinates  should 
form  a  roughly  rectangular  mesh  (fig.  4).  The  method  uses  a  linear  interpola¬ 
tion  scheme  (ref.  6)  to  determine  an  approximate  ray-surface  intersection 
point.  Once  a  "near  solution"  point  is  found,  the  method  calls  the  spline  fit 
routine  in  the  two  orthgonal  directions  (Y  and  Z  in  fig.  6)  to  calculate  a 
more  accurate  interception  point.  The  subroutine  checks  the  interception 
points  and  the  changes  in  slope  using  the  spline  fit  in  both  orthogonal  direc¬ 
tions.  Should  they  be  outside  the  input  tolerances,  an  error  message  is 
printed.  However,  the  calculation  continues. 

The  surface  normal  is  calculated  in  a  manner  similar  to  that  of  the 
general  axi symmetrl c  case.  The  differences  are  that  the  outputs  from  the 
cubic  spline  fits  include: 

dX/dY  and  dX/dZ 

Thus 

tan  £  =  1 /(dX/dY)  (21) 


and 


tan  K  =  1 /(dX/dZ)  (22) 


Now 


tan  a  =  c/d  (23) 


and 


c  =  (a2  +  b2)1/2  (24) 


where 

a  =  b  tan  £  (25) 


9 


and 


d  =  b  tan  < 


where  a,  b,  c,  and  d  are  distances  defined  in  figure  6,  so 


tan  a  =  b(tan2£  +  l)1/2/(b  tan  <) 


tan  a  =  (tan2E  +  1 ) 1  /2 / tan  < 


Once  the  value  of  a  can  be  determined,  v  and  £  can  be  evaluated  in  a  man¬ 
ner  similar  to  the  general  axi symmetric  case  where: 

cos  v  =  1/(1  +  tan2o  +  tan2v)^2  (29; 


cos  C  *  (1  -  cos' 


cos2v) 1 /2 


It  is  possible  to  describe  the  interception  point  and  the  surface  normals 
at  that  point.  Now  the  changes  in  probe  volume  geometry  can  be  calculated. 

TOLERANCES 

The  computer  program  requests  user  input  of  four  tolerances.  There  are 
two  types  of  tolerances.  The  first  type  specifies  how  near  parallel  to  a  coor¬ 
dinate  axis  a  direction  vector  need  be  to  be  considered  parallel  to  that  axis. 
This  simplifies  some  calculations  and  can  improve  accuracy  by  eliminating  divi¬ 
sions  by  numbers  near  zero.  The  other  type  of  tolerance  is  used  to  specify 
the  solution  accuracy  of  the  window  surfaces  during  the  calculations  of  the 
general  smooth  windows. 

All  requested  tolerances  have  default  values. 

The  first  tolerance  requested,  ' TL 1 ' ,  is  used  during  the  calculation  of 
the  direction  cosines  of  the  nonrefracted  laser  beams.  In  later  cal cuiations , 
divisions  by  the  direction  cosines  of  the  laser  beams  are  required.  The  toler¬ 
ance  *  TL 1  *  specifies  at  what  point  the  beam  is  considered  to  be  parallel  to  a 
coordinate  axis.  This  insures  that  division  by  a  very  small  number  during  the 
calculation  of  the  beam  surface  intercept  point  does  not  occur. 


Tolerance  ’ TL2 '  restricts  the  potential  error  generated  in  the  solution 
of  the  general  surface  outward  normal.  Small  errors  in  these  values  can  make 
noticeable  differences  in  the  final  solution.  To  limit  the  potential  error, 
the  slopes  of  the  four  points  bracketting  the  proposed  interception  point  are 
calculated.  The  changes  in  the  slopes  in  the  two  coordinate  directions  are 
calculated  and  these  changes  are  required  to  be  less  than  'TL2.'  This  effec¬ 
tively  limits  the  maximum  possible  error  in  the  surface  normal  of  the  inter¬ 
cept.  If  this  convergence  criteria  is  not  met,  a  finer  grid  describing  the 
window  surface  is  generated  and  the  calculation  is  repeated. 

The  third  tolerance,  'TL3',  restricts  the  possible  error  of  the  intercept 
point  of  the  incoming  light  beam  and  the  window  surface.  The  intercept  point 
calculated  using  a  linear  interpolation  between  the  points  bracketting  the 
intercept  region  is  compared  with  the  intercept  point  calculated  using  a  cubic 
spline  interpolation.  If  the  difference  is  within  the  tolerance  ’ TL3 1  the  pro¬ 
gram  continues.  If  this  criteria  is  not  met,  a  finer  grid  describing  the  sur¬ 
face  is  generated  and  the  intercept  calculation  is  repeated. 

The  last  input  tolerance,  1 TL4 1 ,  is  similar  to  the  first  since  it  speci¬ 
fies  whether  or  not  a  direction  cosine  is  near  parallel  to  the  X-axis  or  lies 
in  the  X  to  Y  or  X  to  Z  planes.  Those  special  cases  allow  simplified 
calculations  and  use  of  the  simplified  calculations  insures  that  division  by  a 
very  small  number  does  not  occur. 

All  default  tolerances  result  in  sufficiently  accurate  solutions  for  the 
types  of  window  in  use  at  the  Lewis  Research  Center.  For  comparison,  note 
that  the  typical  window  at  Lewis  Research  Center  might  have  a  change  in  slope 
of  2.3°  per  centimeter  in  a  direction  perpendicular  to  the  incoming  Laser  Beam. 


EXAMPLES 

A  number  of  test  cases  for  the  flat  plate  and  cylindrical  window  were  run 
and  the  results  are  available  in  reference  8.  A  comparison  of  output  with  an 
unpublished  algorithm  and  also  with  hand  calculations  show  an  exact  agreement 
with  both  the  flat  plate  and  simple  cylinder  test  cases. 

Since  the  window  shapes  for  the  general  axisymmetric  and  smooth  general 
windows  can  only  be  approx imated ,  ideal  comparisons  are  difficult.  In  results 
not  given  here,  the  use  of  the  general  axisymmetric  technique  to  calculate  win¬ 
dow  effects  of  a  cylindrical  window  results  in  an  identical  solution  with  the 
cylindrical  window  solution  and  using  the  smooth  general  technique  to  calcu¬ 
late  the  effects  of  a  flat  plate  window  also  results  in  an  exact  comparison 
with  the  flat  plate  solution. 

The  general  axisymmetric  case  was  also  run  to  calculate  the  window  effects 
of  an  actual  compressor  window.  Figure  7  shows  the  input  geometry.  The 
bisecting  line  is  moved  out  of  the  X  to  Z  and  X  to  Y  planes.  This  gen¬ 
eral  axisymmetric  case  was  also  compared  with  two  outputs  run  using  the  general 
smooth  window  algorithm.  These  two  general  smooth  cases  varied  only  with 
respect  to  the  orientation  of  the  window  and  lens  geometry  to  the  coordinate 
system.  The  first  orientation  run  was  similar  to  the  general  axisymmetric 
orientation  shown  in  figure  7.  The  second  orientation  is  shown  in  figure  8 
and  simply  results  from  rotating  the  problem  geometry  about  the  Y-axis.  The 


window  geometry  is  given  in  Table  1(a)  and  the  other  relevant  parameters  are 
given  in  Table  1(b) . 

A  comparison  of  outputs,  provided  in  figures  9  to  11,  showed  very  good 
agreement  between  the  three  cases. 

Figure  9  shows  the  uncrossing  results,  in  fringes,  versus  6  for  the 
three  cases.  Maximum  difference  between  cases  is  slightly  over  three  fringes 
at  S  equal  90°.  This  is  a  relatively  large  number  of  fringes;  however,  it 
occurs  in  a  geometry  where  the  total  uncrossing  is  a  minimum.  Over  most  of 
the  range  of  comparison,  the  difference  is  less  than  two  fringes.  From  a  prac¬ 
tical  standpoint,  the  agreement  for  all  three  cases  is  excellent. 

Figure  10  shows  the  change  in  beam  crossing  angled  versus  6  for  the 
three  test  cases.  Once  again,  the  agreement  is  very  good  with  the  maximum  dif¬ 
ference  in  the  three  conf igurations  less  than  0.017°.  This  maximum  difference 
results  in  velocity  errors  of  less  than  0.3  meters  per  second.  Note,  however, 
that  the  maximum  calculated  velocity  error  due  to  the  change  in  crossing  angle 
from  refraction  effects  is  calculated  to  be  about  1  meter  per  second  for  these 
test  cases. 

Figure  11  shows  the  change  in  probe  volume  position  for  each  of  the  three 
cases.  The  figure  shows  the  relative  movement  in  each  coordinate  direction 
and  the  total  relative  movement.  Movement  is  defined  as  the  change  in  probe 
volume  position  from  the  position  it  would  occupy  if  there  were  no  intervening 
window.  Relative  movement  is  that  distance  divided  by  the  window  thickness. 
Agreement  is  excellent  for  total  movement  and  X  direction  movement  with  a 
maximum  difference  being  less  than  10  percent  or,  in  absolute  differences, 
0.00635  cm  (window  thickness  of  0.3175  cm).  The  Y  and  Z  direction  move¬ 
ments  show  excellent  agreement  between  the  two  test  cases  with  the  same  win¬ 
dow/coordinate  system  orientation.  However,  when  the  window  surface  is 
rotated  to  an  orientation  that  is  not  approximately  in  the  Y  to  Z  plane, 
the  agreement  is  not  as  good.  Nevertheless,  for  both  directions,  the  total 
difference  in  movement  between  test  cases  is  only  about  0.00225x  window  thick¬ 
ness  or  about  7.15xlO-^  cm  for  a  0.3175  cm  thick  window.  This  is  equivalent 
to  less  than  one  third  of  the  probe  volume  diameter. 

This  set  of  test  cases  demonstrates  two  important  points.  The  first  is 
to  use  the  surface  description  that  most  accurately  reflects  the  actual  window 
geometry.  If  the  window  can  be  described  using  a  general  axisymmetric  surface, 
this  is  the  most  desirable  method  of  input.  If  the  window  can  be  specified  as 
a  flat  plate  or  simple  cylinder,  these  are  the  most  desirable  datasets.  The 
second  point  is  that,  if  the  surface  to  be  considered  is  best  described  as  a 
smooth  general  surface,  the  orientation  of  the  window  surface  normal  should  be 
kept  as  nearly  parallel  to  the  X-axis  as  possible. 


CONCLUSIONS 

A  FORTRAN  computer  program  has  been  developed  to  solve  for  the  movement, 
beam  uncrossing,  and  bisecting  angle  change  of  a  probe  volume  of  a  laser  fringe 
anemometer  when  it  must  traverse  an  intervening  region  of  different  refractive 
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APPENDIX  A 


INPUTS  AND  OPERATION  OF  THE  PROGRAMS 

As  currently  configured,  two  separate  computer  codes  are  necessary  for  a 
complete  solution.  Both  programs  are  interactive  and  require  user  inputs. 

Both  programs  use  FORTRAN  coding  for  the  terminal  interactive  operations  and 
for  manipulating  the  necessary  data  files. 

The  first  program  sets  up  an  input  data  set  to  be  run.  This  allows  a  per 
manent  file  of  different  input  data  sets  to  be  maintained.  The  input  data 
file  structure  is  shown  in  figure  A1 .  A  glossary  of  terms  is  provided  in 
Appendix  B. 

Lines  4  through  6  are  optional.  Should  nop2  be  set  to  zero,  the  program 
will  generate  output  for  beam  plane  orientations  every  10°  from  0  to  170.  In 
this  instance,  no  entries  need  to  be  made  for  lines  4  through  6  and  the  user 
will  be  immediately  prompted  for  line  7. 

Starting  with  line  9,  the  input  varies  depending  on  the  nature  of  the  win 
dow  to  be  described.  For  a  general  window,  the  "upper"  or  outer  surface  varia 
bles  are  read  in  one  at  a  time  increasing  in  the  * Y '  or  leading  index  first, 
as  shown  in  figure  4.  When  finished  with  the  window  upper  surface,  continue 
with  the  lower  or  inner  surface. 

Inputs  for  the  flat  plate  are  done  in  the  same  manner  except  that  only 
four  points  per  surface  (total  of  eight  for  the  window)  should  be  entered. 

For  the  cylindrical  window,  the  only  input  past  card  8  is  the  outer  sur¬ 
face  radius. 

For  the  general  axisymmetric  window,  type  in  the  number  of  Z  to  R 
points  to  describe  the  upper  surface  and  the  Z  to  R  values. 

Figure  A2  shows  the  coordinate  system  and  variables  used  to  specify  the 
problem  geometry.  Angles  should  be  input  as  degrees  and  distances  can  be  spec 
ified  in  any  consistent  set  of  units. 

The  second,  or  beam  uncrossing  program  requires  the  file  name  containing 
the  input  data. 

The  beam  uncrossing  program  then  queries  the  operator  to  insure  that  the 
input  parameters  are  correct.  This  also  allows  modification  of  all  inputs 
other  than  the  basic  window  geometry. 
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All  Integer 
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INPUT  DICTIONARY 


fields  are  15  format  and  all  real  fields  are  F10  format. 

an  array  of  angles  describing  the  angle  between  the 
plane  the  two  incoming  beams  lie  in  and  a  plane  contain¬ 
ing  the  lens  center-z  axis  (axis  of  rotation),  see  fig¬ 
ure  A2b,  In  degrees. 

specifies  the  magnitude  of  the  second  derivative  at  the 
end  points  of  the  cubic  spline  fit  when  the  window  to 
be  modeled  is  a  general  axisymmetric  or  general  smooth 
window.  (default=.5) 

the  bisecting  angle  of  the  beams  in  degrees 

number  of  points  in  the  "y"  direction  to  describe  sur¬ 
faces  1,  2;  01  for  cylindrical,  02  for  flat  plate,  and 
10(max)  for  general  surface 

number  of  points  in  the  "z"  direction  to  describe  sur¬ 
faces  1,  2;  02  for  flat  plate,  10(max)  for  general 
surface 

index  of  refaction  of  the  medium  before  and  after  trav¬ 
ersing  the  window;  (default  =  1.00029,  air) 

number  of  Z  to  R  coordinate  pairs  used  to  describe 
a  general  axisymmetric  window 

index  of  refraction  of  the  window;  (default  =  1.6,  nomi¬ 
nal  glass) 

flag  to  define  the  type  of  window  to  be  input;  0  -  gen¬ 
eral  surface,  1  -  flat  plate,  2  -  cylindrical,  and  3  - 
general  axisymmetric 

defines  whether  the  series  of  angles  (delta)  is 
defaulted  or  specified;  0  -  defaulted,  gives  delta 
increments  of  10°  from  0  to  170,  1  to  18  indicates  nop2 
number  of  delta  angles  to  be  input 

outer  surface  radius,  general  axisymmetric  or  cylindri¬ 
cal  surface. 

distance  from  the  center  of  the  focusing  lens  to  the 
axis  of  rotation  (z  axis) 

distance  of  the  probe  volume  (lens  focal  point  assuming 
no  intervening  window)  to  the  axis  of  rotation  (z  axis) 
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thk  window  normal  thickness;  not  used  for  general  surface 

Input. 

thl  angle  between  the  x  axis  and  a  perpendicular  line  con¬ 

necting  the  z  axis  and  the  center  of  the  focusing  lens 

thpv  angle  between  the  x  axis  and  a  perpendicular  line  con¬ 

necting  the  z  axis  and  the  probe  volume  (focal  point  of 
the  lens) 

til  tolerance  used  in  subroutine  'ppos.'  Defines  whether 

or  not  the  plane  in  which  the  two  incoming  laser  beams 
lie  is  parallel  to  any  coordinate  plane.  This  speci¬ 
fies  the  nature  of  the  calculation  of  the  beam/focusing 
lens  plane  intercept  calculation,  (default  is  .0017rad) 

t 1 2  tolerance  in  subroutine  1 sf i n ' .  Determines  a  success¬ 

ful  solution  to  the  general  surface/laser  beam  inter¬ 
cept.  This  criteria  restricts  the  maximum  error 
between  the  actual  surface  normal  direction  and  the  cal¬ 
culated  value.  Default  is  set  to  cos  (angle  error)c.01. 

t 1 3  tolerance  in  subroutine  ' sf i n ‘ .  Determines  a  success¬ 

ful  solution  to  the  general  surface/laser  beam  inter¬ 
cept.  This  is  a  limit  between  the  difference  in  calcu¬ 
lated  intercept  points  using  the  'linint'  (linear  inter¬ 
polation)  approximation  of  the  window  surface  and  the 
'splint'  (spline  interpolation)  approximation  of  the 
window  surface.  It  is  assumed  that  when  the  difference 
is  less  than  < 1 1 3* th! .) ,  the  surface  description  is  accu¬ 
rate  for  an  acceptable  answer.  Default  is  0.0015. 

tl4  tolerance  in  subroutine  'sfin'.  Determines  if  any  sur¬ 

face  direction  vector  is  near  parallel  to  the  x-axis 
or  liens  in  the  X  to  Y  or  X  to  Z  planes.  Maxi¬ 
mum  off  parallel  or  out-of-plane  error  for  the  direc¬ 
tion  cosine  is  t!4.  Default  is  .0017. 

xl ,yl ,zl(1-iyl ,1-izl )  x,  y,  z  coordinates  describing  the  upper  surface  for 

the  flat  plate  or  general  window  case. 

x2,y2,z2(l-1y2,l-1z2)  x,  y,  z  coordinates  describing  the  lower  surface  for 

the  flat  plate  or  general  surface  case 

wvl  Laser  beam  wavelength  in  Angstroms.  Used  for  calculat¬ 

ing  the  uncrossing  in  fringes,  (defaul t=51 45A) 


zcl ( 1 -nr) 


z  coordinates  (along  the  axis  of  symmetry)  correspond¬ 
ing  to  rcl  for  the  general  axi symmetri c  case 

axial  distance  of  focusing  lens  center 

axial  distance  of  probe  volume  (focal  point) 
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SYMBOLS 

f  prove  volume  fringe  crossing  frequency 

n  index  of  refraction 

N  outward  surface  normal  at  point  of  laser  beam/surface  intersection 
V  fluid  (particle)  velocity  vector 
Vi  laser  beam  direction  vector,  i  =  1  or  2 
o,3, y  direction  angle  for  incoming  laser  beam 

v.C.a  direction  angles  for  outward  surface  normal  vector  (fig.  5  and  6) 

X  1 ight  wavelength 

K  angle  between  a  surface  tangent  line  in  the  X  to  Z  plane  and  the  Z 

axis  (fig.  5  and  6) 

E  angle  between  a  surface  tangent  line  in  the  X  to  Y  plane  and  the 

Y  axis  (fig.  5  and  6) 

e  beam  crossing  angle  in  probe  volume 

4>  angle  between  the  outward  surface  normal  and  the  incoming  laser  beam 


Superscript 

1  indicates  the  region  after  the  laser  beam  has  crossed  the  surface 

boundary 

"  indicates  a  coordinate  system  centered  on  the  surface/ray  intercept 
point.  Used  in  figures  4  and  6. 
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SAMPLE  INPUTS  AND  OUTPUTS 

Figures  01  through  D4  show  the  sample  Inpucs  for  flat  plate,  simple  cylin¬ 
der,  general  axlsymmetrlc,  and  general  smooth  windows,  respectively.  In  these 
figures,  all  user  Inputs  are  underlined  and  an  asterisk  Is  found  at  the  right 
hand  margin  of  these  lines. 

Notice  that  the  data  set  up  program  allows  the  option  to  change  any  input 
data  line.  The  line  number  of  the  line  that  is  to  be  changed  is  typed  in, 
then  the  entire  line  Is  retyped. 

Figure  D5  shows  a  sample  output  from  the  main  optics  code.  This  output 
Is  self  explanatory. 
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WARNING  MESSAGES 

A  number  of  warning  messages  may  be  printed  during  the  course  of  a  calcu¬ 
lation.  They  are,  for  the  most  part,  only  advisory  in  nature  with  the  calcula 
tion  continuing.  However,  they  may  indicate  an  error  in  the  input  data.  The 
possible  warning  messages  are: 

1.  "WARNING:  FROM  RTRACE,  INTERNAL  REFLECTION  OCCURS  IN  THIS  CONFIGURATION." 

It  is  possible  for  an  optics  window  conf iguration  to  occur  where  the  beam  leav 
ing  the  window  intercepts  the  window  surface  at  an  angle  greater  than  Bragg's 
angle.  Should  this  occur,  an  internal  reflection  of  the  beam  will  occur.  The 
calculation  stops  and  output  from  this  configuration  should  be  considered 
Incorrect.  This  situation  should  only  occur  with  general  smooth  and  general 
ax i symmetric  windows. 

2.  "WARNING:  FROM  LINT.  LINES  DO  NOT  CROSS  IN  SPACE." 

This  warning  is  printed  out  when  calculations  indicate  that  the  two  incoming 
beams  lie  in  the  same  plane  but  intercept  point  calculations  do  not  indicate 
the  beams  cross.  This  message  indicates  an  input  data  error. 

3.  "WARNING:  FROM  CYLINT,  LINE  DOES  NOT  INTERCEPT  CYLINDER  WALL." 

This  message  indicates  an  input  geometry  where  the  beams  do  not  intercept  the 
cylinder  window  outer  surface.  This  message  indicates  an  input  data  error. 
Calculation  terminates. 

4.  "WARNING:  FROM  CYLINT,  THE  CALCULATIONS  ON  Z  POSITION  DO  NOT  AGREE." 

This  message  indicates  an  input  data  error  or  severe  geometry  constraints  such 
as  an  intercept  near  the  surface  tangency  point  or  an  excessive  <J>.  Calcula¬ 
tion  continues. 

5.  "WARNING:  FROM  CYLINT,  CALCULATED  INTERCEPT  X  AND  Y  DOES  NOT  LIE  ON  R. 
ALTERNATE  CALCULATION  STARTED." 

This  message  indicates  severe  geometry  constraints  where  the  surface  intercept 
point  lies  near  the  surface  tangency  point.  The  alternate  calculation  should 
provide  greater  accuracy  for  this  geometry.  Calculation  continues. 

6.  "WARNING:  FROM  AXI ,  NO  CONVERGENCE  IN  'AXI'." 

The  calculation  of  the  intercept  point  of  the  general  axi symmetr i c  surface  and 
an  incoming  ray  is  an  iterative  procedure.  Should  the  difference  between  the 
calculated  line  intercept  and  the  surface  intercept  point  be  greater  than  TLl 
after  50  iterations,  this  message  is  printed.  Calculation  continues;  howevc  . 
the  output  data  may  be  in  error.  This  message  indicates  an  input  data  error 
or  severe  surface  curvature  problems. 


7.  "WARNING:  FROM  SURLNT,  SURLNT  FAILS." 

This  message  Indicates  that  the  Intercept  of  the  incoming  light  ray  with  the 
linear  approximation  to  the  general  surface  fails  t->  converge  within  76  itera 
tions.  Calculation  continues;  however,  a  possible  problem  may  exist  with  the 
input  surface  data. 

8.  "WARNING:  FROM  SFIN,  CHANGES  IN  THE  SURFACE  OR  THE  SURFACE  NORMAL  ARE 
RAPID,  SFIN  DOES  NOT  CONVERGE.  SOLUTION  ACCURACY  IS  NOT  GUARANTEED." 

This  message  indicates  a  severe  window  geometry  resulting  in  a  solution  which 
may  be  in  error.  Calculation  continues. 

9.  "WARNING:  FROM  FINEE,  MORE  THAN  ONE  DIRECTION  COSINE  IS  EQUAL  TO  ONE." 

This  message  should  not  occur  and  indicates  an  input  data  error  that  is  proba 
bly  In  the  lens-PV  geometry.  Calculation  terminates. 

10.  "WARNING:  FROM  FINEE,  THE  BEAMS  ARE  PARALLEL  AND  DO  NOT  CROSS." 

This  message  indicates  that  the  window  geometry  has  resulted  in  two  parallel 
beams.  There  will  be  no  probe  volume  and  no  minimum  crossing  distance.  Calc 
latlon  terminates. 

11.  "WARNING:  FROM  FINEE,  DMIN='  ',  DMN= 1 

The  minimum  crossing  distance  is  calculated  in  two  ways.  Should  the  results 
be  different,  this  advisory  only  message  is  printed.  A  large  difference  in 
values  would  indicate  an  input  data  error. 
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